Controlled collisions of a single atom and ion guided by movable trapping potentials 
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We consider a system composed of a trapped atom and a trapped ion. The ion charge induces in 
the atom an electric dipole moment, which attracts it with an r -4 dependence at large distances. 
In the regime considered here, the characteristic range of the atom-ion interaction is comparable or 
larger than the characteristic size of the trapping potential, which excludes the application of the 
contact pseudopotential. The short-range part of the interaction is described in the framework of 
quantum-defect theory, by introducing some short-range parameters, which can be related to the 
s-wave scattering length. When the separation between traps is changed we observe trap-induced 
shape resonances between molecular bound states and vibrational states of the external trapping 
potential. Our analysis is extended to quasi-one-dimensional geometries, when the scattering exhibit 
confinement-induced resonances, similar to the ones studied before for short-range interactions. For 
quasi-one-dimensional systems we investigate the effects of coupling between the center of mass and 
relative motion, which occurs for different trapping frequencies of atom and ion traps. Finally, we 
show how the two types of resonances can be employed for quantum state control and spectroscopy 
of atom-ion molecules. 



PACS numbers: 32.80.Pj, 34.90.+q 
I. INTRODUCTION 



Techniques developed in atomic physics during the last 
two decades allow the preparation of single atoms and 
ions in the laboratory. Single neutral atoms can be stored 
and laser cooled in far off-resonance laser traps (FORT) 
0, 0, and arrays of atoms can be prepared in an op- 
tical lattice via a Mott insulator phase, where exactly 
one atom is stored per lattice site [3[ . Similar techniques 
are being pursued in the context of atom chips [|| . Fur- 
thermore, single ions and arrays of ions can be stored 
in Paul and Penning traps, and sideband laser cooling 
allows to prepare ions in the vibrational ground state 
of these trapping potentials [f|. The controlled prepa- 
ration, manipulation and measurement of electronic in- 
ternal and motional states of laser driven single atoms 
and ions, combined with the possibility of controlled en- 
tanglement of atoms and ions provides the basic ingredi- 
ents to for investigations of fundamental aspects of quan- 
tum mechanics, and applications such as high-precision 
measurements and quantum information. In quantum 
information processing entanglement of qubits stored in 
internal states of single atoms or ions is achieved by con- 
trolled interactions between particles, either in the form 
of switchable qubit-dependent two-particle interactions 
or via an auxiliary collective mode of the systems, which 
serves as a quantum data bus. In particular, for neutral 
atoms cold controlled collisions of atoms stored in mov- 



able spin-(i.e. qubit-) dependent optical lattices has been 
proposed as a means to entangle atomic pairs, and has 
experimentally implemented to generate iV-atom cluster 
states Q. 

In the present work we study the controlled cold colli- 
sion of a single atom and a single ion, where the atom 
(ion) is prepared in a given motional state of an atom 
(ion) trap, and we move the traps to guide the atom 
and ion wave packet to "collide" for a given time (see 
Fig. [I]). Our focus is the development of a quantum- 
defect formalism for trapped atoms and ions, in a form 
which is convenient for future applications and exten- 
sions, in particular in the context of quantum informa- 
tion processing to swap qubits stored in atoms and ions, 
and the entanglement of these qubits in a controlled colli- 
sion. The atom-ion collision is governed by the potential 
V(r) — > — ae 2 /(2r 4 ) (r — > oo), where a is the dipolar 
polarizability - to be contrasted to a van der Waals po- 
tential V(r) — > — Cs/(r 6 ) which represents the collisional 
interactions between e.g. alkali atoms in their electronic 
ground state. Below we will calculate the dynamics of the 
interacting atom-ion system for a given time-dependence 
of the motion of the trap. This includes the transition to 
excited trap states after the collision, the possible forma- 
tion of an excited atom-ion (molecule) complex, and the 
description of trap induced resonances. 

To study the controlled atom-ion collisions we first 
derive an effective Hamiltonian, describing effective net 
forces acting on the particles moving in rapidly changing 
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FIG. 1: (Color online) Schematic drawing of a controlled 
collision between a single atom and a single ion, whose center- 
of-mass wavepackets are guided by time-dependent atom trap 
and ion trap, respectively. Labels distinguish different phases 
of the process: a) initially particles prepared in the motional 
ground state; b) collision (overlap of the wave packets); c) 
excitation of the motional states in the traps; d) particles in 
some excited states after the collision. 



laser and rf fields. In our approach, we include in the 
Hamiltonian only the asymptotic part of the atom-ion 
potential, whereas the short-range interactions are taken 
into account by imposing appropriate boundary condi- 
tions on the wave function at r — > 0. In the regime of cold 
collisions considered here, the boundary conditions can 
be expressed in terms of a single quantum-defect param- 
eter, independent of the collisional energy and angular 
momentum. The quantum-defect approach is dictated 
by the relatively long-range character of the atom-ion in- 
teraction, which exceeds the typical size of the trapping 
potentials, and, in contrast to the neutral atoms, excludes 
the applicability of the pscudopotcntial. 

Our method to describe the dynamics of the controlled 
collision is based on the application of the correlation di- 
agrams, i.e. energy spectra as a function of the trap 
separation. Such correlation diagrams, widely used in 
quantum chemistry to characterize reactions of diatomic 
molecules, connect in our case the asymptotic vibrational 
states with the molecular and vibrational states at zero 
trap separations. At the intermediate distances the en- 
ergy curves exhibit avoided or diabatic crossings, depend- 
ing on the symmetry of eigenstates and of the coupling 
term in the Hamiltonian. In the atom-ion system the 
avoided crossing can be attributed to the resonances be- 
tween molecular and vibrational states, that appear when 
the energy of a vibrational level, coincide with the energy 
of a molecular state shifted by the external trapping. The 
dynamics in the vicinity of such avoided crossings can be 
accurately described in the framework of the Landau- 
Zener theory. In this paper we determine the level split- 
ting at the avoided crossing, at different trap separations, 
and for different symmetries of the molecular states. In 
this way we can characterize the time scale appropriate 
for adiabatic or diabatic traversing of a given avoided 
crossing. 

In our paper we consider two different geometries of 
the trapping potentials: spherically symmetric traps and 



very elongated cigar-shape traps (quasi-lD traps). The 
former case requires full three-dimensional treatment, 
while the latter one can be described using an effec- 
tive one-dimensional Hamiltonian. In the latter case, 
the tight transverse confinement effectively renormalizcs 
the one-dimensional quantum-defect parameters. In ad- 
dition, very elongated cigar-shape traps exhibit another 
type of resonances, that can be observed at zero-trap sep- 
arations, but for changing ratio of the transverse confine- 
ment to the s-wave scattering length of the atom-ion in- 
teraction. Such confinement-induced resonances appear 
when the energy of the colliding particles coincide with 
energy of a bound state lifted up by the tight transverse 
potential. 

As already outlined above, the counterpart to the idea 
of ultracold controlled trap-guided collisions has been im- 
plemented in optical lattices, where application of spin- 
dependent potentials has allowed to entangle atoms be- 
tween neighboring sites Moreover, resonance phe- 
nomena in trapped systems, similar to the ones consid- 
ered here, have been already thoroughly investigated for 
interactions between neutral atoms. In quasi-lD systems, 
collisions between neutral atoms exhibit confinement- 
induced resonances (Io| . whereas displacement of the 
trapping potentials leads to trap-induced resonances pd| . 

The paper is organized as follows. In section |TT] we 
present our model that we use to describe the controlled 
atom-ion collisions. In particular, in section III Al we in- 
troduce an effective Hamiltonian of a single trapped atom 
and ion. Section Hi Bl discuss the quantum-defect treat- 
ment of the short-range part of the interaction potential. 
In sections fll Cl andl ll Dl we discuss our model for two par- 
ticular geometries of the trapping potentials, spherically 
symmetric traps of the same frequency and elongated 
cigar-shape traps with the same transverse trapping fre- 
quency, respectively. We extend the quantum-defect 
treatment to quasi-lD systems in section Til El Our ap- 
proach to describe the dynamics of the atom-ion system is 
discussed in section Hi Fl Section Hill presents the results 
of our calculations. We start from discussing the proper- 
ties of the scattering in quasi-lD systems in the presence 
of long-range r~ A potential, in particular, the problem 
of determining of ID quantum-defect parameters. Sec- 
tions IIII Bl and IIII CI present adiabatic energy curves and 
adiabatic eigenstates in, respectively, ID and 3D traps 
with the same trapping frequencies for the atom and the 
ion. The dynamics in the vicinity of the avoided crossing 
is discussed in section llll Dl where we apply the Landau- 
Zener theory and semiclassical methods to address this 
problem. Section IIII El discuss most complicated case of 
different trapping frequencies, focusing on the ID atom- 
ion motion. The outlook and final conclusions are pre- 
sented in section IIVI Finally appendix [X] presents the 
microscopic derivation of the effective Hamiltonian and 
appendix [B] describes the details of the pseudopotential 
treatment of the scattering in quasi-lD traps. 
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II. BASIC SETUP AND MODEL 

We consider a system consisting of a single atom and 
single ion, stored in their respective trapping potentials. 
Single atoms can be trapped in experiments with optical 
tweezers [l], 0] or optical lattices Q . Such potentials are 
created in far-detuned laser fields due to the AC Stark 
effect. In another experimental technique, on atom chips 
H| , single atoms can be trapped in microtraps created by 
electric and magnetic fields around wires and electrodes. 
On the other hand, single ions can be confined in radio- 
frequency traps, which use a rapidly oscillating electric 
field to create an adiabatic trapping potential for the ion 
charge 0. Typically, the trapping potentials close to 
the trap center are with a good approximation harmonic, 
and in our approach we consider the simple picture of 
the atom and ion confined in the harmonic potentials 
with frequencies u a and u>i, respectively, with d denoting 
the displacement between the traps (c.f. Fig. [2Ja). We 
assume that the atom and ion are prepared initially in a 
given trap state, e.g. the ground state of the potential, 
and that the atom and ion are brought to collision by 
overlapping their center-of-mass (COM) wave packets, 
i.e. we change d(t) as a function of time as outlined in 
Fig. [T] 

In studying this collision with COM wave packets 
guided by movable trapping potentials we will study two 
cases of (i) three-dimensional (3D) dynamics, where both 
traps are assumed to be spherically symmetric, and (ii) 
the case of quasi-lD traps, i.e. elongated, cigar-shaped 
traps which can be described as an effectively ID colli- 
sion dynamics. This second case is conceptually simpler, 
but also provides collisional features specific to ID sit- 
uation (e.g. the appearance of ID confinement-induced 
resonances, reminiscent of resonances in [Io|]). The val- 
ues of trapping frequencies in atom and ion traps are 
typically quite different, thus it is natural to consider 
Lj a ^ u>i. For such condition, however, the COM and 
relative motion are coupled, which significantly compli- 
cates the theoretical description. Apart from this general 
regime, we consider also the special case of u a — Wj, when 
COM and relative motions can be separated. In this case 
we perform the full diagonalization for the 3D problem, 
which is numerically very demanding in the general case 

Of lQ a ^ 10 i. 

At large distances, the interaction between atom and 
ion is given by a potential V(r) ~ —ae 2 / (2r 4 ), where a is 
the dipolar polarizability. This potential originates from 
the attraction between the ion charge and the electric 
dipole that it induces on the atom [24| ■ At short distances 
the interaction is more complicated; however, it turns out 
that a detailed knowledge of the form of core region of the 
potential is not necessary, and that it can be subsumed in 
our model by introducing a set of (energy independent) 
quantum-defect parameters. 

Let us focus now on characteristic scales that can be 
associated with the atom-ion interaction. For the sake 
of clarity table Q] summarizes all the length scales that 
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FIG. 2: (Color online) Panel (a): Schematic drawing of the 
considered setup consisting of a trapped atom and ion. Panels 
(b) and (c): Long and short-range regimes, respectively, in the 
controlled collisions between a trapped atom and ion. 



are introduced in the course of the paper. We define 
its characteristic length as R* = y ae 2 /j,/H 2 , and char- 
acteristic energy as E* = h 2 / (2/i(R*) 2 ) , where /i is the 
reduced mass. Other characteristic lengths and energies 
are given by the ground states of the trapping poten- 
tials. We associate to the trapping frequencies the har- 
monic oscillator lengths for atom (ion) : l a = y/h/ (m a cj a ) 

(li = \Jh/(miUJi)). For optical and rf traps, the trapping 
frequencies are typically uj a = 2ir x fO — fOOkHz and 
u>i = 2ir x 0.1 — fOMHz. In table HT1 we collected some 
example values of R*, U, l a , and E* for a few systems 
of alkali atoms and alkaline earth ions. One can observe 
that the range of the polarization potential is compara- 
ble or larger than the size of the trapped ground state. 
Therefore, for the trapped particles the polarization in- 
teraction cannot be replaced with the standard contact 
pscudopotcntial familiar from s-wavc atom-atom scatter- 
ing in the context, e.g., of Bose-Einstein condensation. 

In the interactions between trapped atom and ion, one 
can distinguish two different regimes depending on the 
relative distance d between traps: (i) d ^> Ri 1 R ai (ii) 

d < Ri ~ R a , where R v = (ae 2 /m„w^) 1 6 for v = i,a 
is some characteristic distance at which the atom-ion in- 
teraction becomes comparable to the trapping potential 
of the ion and atom, respectively. Because of the weak 
dependence of R v on ui v and m v , the characteristic dis- 
tances Ri and R a are roughly the same. The two regimes 
(i) and (ii) correspond to the cases where the motion in 
the traps is weakly or strongly affected by the atom-ion 
interactions, respectively The discussed regimes are il- 
lustrated schematically in panels (b) and (c) of Fig. [2l 
and they correspond to snap shots in the collision pro- 
cess when the wavepackets of the particles do not over- 
lap / do overlap during the controlled collision. In the 
regime (i) of large distances, the atom-ion interaction can 
be treated pcrturbatively as a distortion of the trapping 
potential. In this limit the system can be described in 
terms of two coupled harmonic oscillators, and all the dy- 
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Definition 



Description 



Effective Hamiltonian 



R* = y/ae 2 fi/h 2 



la = y/h/(m a Ula) 
I = y/h/(fjw) 



l± = y/h/(llLJj_) 



2 \ 1/6 



Ri 



Ra — 



Rtc\ — 



2 \ 1/6 



2 \ 1/6 



R 



- ( ae2 



2 \ 1/6 



Rid = max(i?x, l±) 
Ro 



characteristic length of the polarization 
potential 

harmonic oscillator length of ion trap 

harmonic oscillator length of atom trap 

harmonic oscillator length for relative 
degrees of freedom (axial direction in 
case of quasi ID traps) 

harmonic oscillator length in the trans- 
verse direction for relative degrees of 
freedom 

distance at which the polarization po- 
tential becomes equal to the ion trapp- 
ping potential 

distance at which the polarization po- 
tential becomes equal to the atom 
trappping potential 

distance at which the polarization po- 
tential becomes equal to the trappping 
potential for relative degrees of freedom 

distance at which the polarization po- 
tential becomes equal to the trapping 
potential in the transverse direction for 
relative degrees of freedom 
boundary of the quasi- ID regime 
size of the core region of the atom-ion 
complex 



TABLE I: Definitions of the length scales used throughout 
the paper. 



We adopt the following time-dependent Hamiltonian 
to describe the system of a single trapped atom and a 
single trapped ion (c.f. Fig. [5]): 



H(t) 



u—i,a 



2m 



Vl , 1 2/ 

+ -m v u) v {z v 



d v {t)f 



\™> v J{ v pl 



(1) 



Here the label i (a) refers to the ion (atom) respectively, 
p and r are the momentum and position operators, d v {t) 
denotes the positions of the atom and ion traps, respec- 
tively, that can be controlled in the course of dynamics, 
and p 2 = x 2 + z 2 . We assume that the trapping po- 
tentials are axially symmetric and displaced along the 
axis of symmetry The trapping frequencies are denoted 
by lo v and ui_ u for the axial and transverse directions, 
respectively Finally, V(r) denotes the interaction po- 
tential between the atom and the ion. At large distances 
the main contribution to this interaction comes from the 
polarization of the atomic cloud, V(r) ~ — ae 2 /(2r 4 ). 

A microscopic derivation of the Hamiltonian |T]) is pre- 
sented in appendix [A} The basic idea of this derivation 
is based on the application of the Born-Oppcnhcimcr ap- 
proximation to the motion of the outer shell electrons, 
and followed by a time averaging over fast time scale of 
the rf and laser frequencies. In this picture, the Hamilto- 
nian (fTJ) represents effective net interactions felt by parti- 
cles moving in the rapidly changing time-dependent po- 
tentials. The interaction potential V(r) can be identi- 
fied with the adiabatic Born-Oppcnhcimer curve of the 
electronic ground-state, that depends on the difference 
between the atom and ion COM coordinates. 





R*(a ) 


k(ao) 


la (do) 


E*/h (kHz) 


40 Ca+ + 87 Rb 


3989 


300 


644 


4.143 


9 Be + + 87 Rb 


2179 


632 


644 


46.59 


40 Ca+ + 23 Na 


2081 


300 


1252 


28.56 



TABLE II: Characteristic distance R* , characteristic energy 
E* , harmonic oscillator length U for ion trap of Wj = 1MHz, 
and harmonic oscillator length l a for atom trap of ui a = 
100kHz, for some combinations of alkali earth ions and al- 
kali atoms. 



namics can be solved analytically. We discuss this case 
in more details in section IIII Fl On the other hand, in 
the regime (ii) the description is more difficult, since it 
requires inclusion of the short-range part of the interac- 
tion, and full treatment of the long-range r -4 part. In 
the present paper we focus mainly on the latter case. 



B. Quantum-defect theory 

We denote by Ro the characteristic distance at which 
V(r) starts to deviate from the asymptotic r~ 4 law. 
Typically, Rq is the size of the core region of the ion- 
atom complex, and is much smaller than all the other 
length scales in our problem. Therefore, we can describe 
the short-range part of interaction in the spirit of the 
quantum-defect theory. For Rq <C R* , the parameters 
describing the short-range potential become independent 
of energy and angular momentum, since for r <C R* the 
interaction potential is much larger than typical kinetic 
energies and heights of the angular-momentum barrier. 
This feature is a key ingredient of quantum-defect theory, 
where the complicated short-range dynamics can be sum- 
marized in terms of few, energy-independent constants 
(phase shifts or quantum defects). 

To analyze the short-distance behavior of the wave 
functions, we omit for the moment the trapping poten- 
tials in the Hamiltonian ([T]), and consider only the part 
describing the scattering of atom from the ion in free 
space. In this case, the relative and COM degrees of 
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freedom are decoupled, and the relative motion is gov- 
erned by the Hamiltonian H = p 2 /2^i + V(r). We apply 
the partial wave expansion to the relative wave function: 
\l/ r ei(r) = J2i Ri( r )Yim(0,4>) with I denoting the angular 
momentum and Y| m (0, 0) the spherical harmonics. For 
r 3> Rq, we set V(r) = —ae 2 /(2r i ) 1 which allows to 
solve radial Schrodingcr equation in terms of the Math- 
ieu functions of the imaginary argument 0, H, Q • In this 
way we obtain the following short-distance behavior of 
the radial wave functions 

Rl(r,k) ~sin(R*/r + ipi(k)), r^^/R*/k, (2) 

where hk is the relative momentum, and (pi(k) are some 
short-range phases. As one can easily verify, the asymp- 
totic solution ^ fulfills the radial Schrodinger equa- 
tion with the energy H 2 k 2 /(2/j,) and centrifugal barrier 
H 2 l(l+1)/ (2/xr 2 ) terms neglected. The short-range phases 
constitute our quantum-defect parameters. For Ro <C R* 
we can assume that fi{k) are independent of the en- 
ergy and angular momentum: <pi(k) = tp, which reduces 
description of the short-range interaction to the single 
quantum-defect parameter <J). In the calculations we re- 
place V(r) by its asymptotic r~ 4 behavior, effectively 
letting Rq — » and imposing boundary condition stated 
by @ with the short-range phase <j>. 

For k = the solution @ becomes valid at all dis- 
tances. Utilizing the fact that zero-energy solution be- 
haves asymptotically as \&(r) ~ 1 — a s /r (r — > oo), we 
can relate the short-range phase to the s-wave scattering 
length a s 

a s = —R* cot ip (3) 

In this way the knowledge of the scattering length a s 
allows us to calculate the quantum defect parameter ip. 

C. Symmetric 3D trapping potential and identical 
trapping frequencies 

In this section we consider 3D dynamics assuming, 
spherically symmetric trapping potentials w„ = lo^ v for 
v = i, a and identical trapping frequencies for atom and 
ion: cui = cu a = u). For such conditions the relative and 
CM motions are decoupled. 

The COM motion is governed by the Hamiltonian of 
a harmonic oscillator with mass M — nii + m a and fre- 
quency uj, while the Hamiltonian of the relative motion 
is given by 

„„, = £!+ i^ (r _ d) *_gl (4) 

where p and r are the relative-motion momentum and 
position operators, respectively, and d = d a — d;. In the 
case of relative motion we define 

R^^iae 2 /^ 2 ) 11 " = (R*) 1/3 l 2/3 (5) 




FIG. 3: (Color online) Contour plot of the total potential 
U(r) — V(r) + ^(iLu 2 (r — d) 2 for the relative motion in the 
atom-ion system. The figure shows the case of spherically 
symmetric trapping potentials with the same trapping fre- 
quencies, separated by d = (0,0, 2R*), and for R* — 21 

(i = v/sTCH). 

as a characteristic length at which the atom-ion inter- 
action is comparable to the trapping potential for the 
relative degrees of freedom, where I = y/h/Jpuj) denotes 
harmonic oscillator length. 

Fig. [21 illustrates a typical potential for the atom-ion 
relative motion. In addition this figure depicts differ- 
ent lengths scales characteristic of our problem. In the 
presented case R* = 21, and the distance between trap 
centers is d = 2R*. 



D. Quasi-ID trapping potentials 

As a second geometry we consider cigar-shape traps 
with a transverse trapping frequency much larger than 
the axial one: lo^ v ^ lo v (quasi-lD traps) for v = i, a. 
For energies smaller than the excitation energy in the 
transverse direction, the motion in the transverse direc- 
tion is frozen to zero-point oscillations, and the dynamics 
takes place along the weakly confined direction. Never- 
theless, the transverse motion plays an important role at 
short distances, effectively renormalizing the short-range 
phase, as we will show in the next section. In this way 
quasi-lD traps offer the additional possibility of tuning 
the interactions, similarly to the case of neutral atoms 
exhibiting confinement-induced resonances [lol ]. 

Here, we consider only some particular situation when 
uj±i = uj± a is the same for atom and ion. This simpli- 
fies our description of the renormalization effects, since 
in this case the transverse COM and relative motions can 
be separated. We expect that our results are also qual- 
itatively valid in the general case of different transverse 
trapping frequencies. 

To obtain an effective ID Hamiltonian, describing the 
evolution of the wave packets along the z axis, we de- 
compose the total wave function of the atom and ion 
into a series over eigenmodes of the transverse part of 
the Hamiltonian. For total energies E < 3hoj± only 
the transverse ground-state mode contributes to the total 
wave function at large distances. In addition, condition 
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\z% — z a \ 3> R±_ assures that the axial and the transverse 
motions are decoupled and the wave function can be writ- 
ten as a product of the axial and the transverse compo- 
nents. Here, R± = (ae 2 / pjJ 2 ) 1 is some characteris- 
tic distance at which the atom-ion interaction becomes 
comparable to the transverse trapping potential, and the 
condition |zj — z a \ ^> Rj_ describes the regime where 
the transverse oscillation frequency is weakly modified 
by the atom-ion interaction. Hence, we have ^(r,, r a ) — ► 

"4>o(Pi,Pa)^lT3{Zi,Z a ) (\Zi ~ Z a \ -» 00 ) for E < 3fejj_. 

Here, Vo is the ground state of the transverse part of 
the Hamiltonian: ip {p h p a ) = e r^^ m ^ +m ^/ 2h /tt 1 ' 2 , 
and ^id(Zi, z a ) denotes the axial part of the wave func- 
tion. Substituting the decomposition over transverse 
modes into the Schrodinger equation with the Hamilto- 
nian l[Tj). and retaining only the lowest transverse mode 
we obtain an effective ID Hamiltonian, that governs the 
dynamics of ^m(zi,z a ) 



H 



ID 



E 



2 

PI 

2m v 



1 



+ -m v it) u (z v - dv) 



V m {\ Zi 



(0) 



Here Vid(|z|) is the ID interaction potential obtained by 
integrating out the transverse degrees of freedom: 



VldG^j - z a \) 



dpidp a \i/Jo(Pi, Pa)\ 2 V(\ri - r a \) 



(7) 



At sufficiently large distances the effective ID interaction 
has similar power dependence as in 3D 



V, 



ID I 



ae 

"2z*' 



\Z > <_L, 



(8) 



where l± = H/(pu)±). Summarizing, in our ID cal- 
culations wc apply the Hamiltonian ^ with the ap- 
proximation |5]). that arc valid for |z, — z a \ Rid = 
max(i?j_, l±). 

Finally in the case of equal longitudinal trapping fre- 
quencies: u>i = Lu a = ui, the relative and the COM de- 
grees of freedom can be separated, and the dynamics is 
described by the Hamiltonian of the relative motion 



H, 



:11D 



2p 



1 



+ -puj 2 (z - dy 



ae 
2^ 



(9) 



where z = z% — z a and p = pi — p a . In the next Section 
we shall solve the problem in the different regimes de- 
scribed before, and investigate the peculiar phenomena 
that arise from the interplay between the trapping and 
the interaction potentials involving the two particles. 



E. Quantum-defect theory in ID 

The quantum-defect treatment of the short-range in- 
teractions can be extended to the ID dynamics, described 



by the Hamiltonian (JH)). In ID the asymptotic behavior of 
the relative wave function at short distances is governed 
by 



^ cl (z,fc) ~ \z\sm{R*/\z\ + ip e {k)), 2 « y/R*/k, 

(10) 

*° el (z,fc) - zsm(R*/\z\ + (p (k)) , z < y/R*/k, 

(11) 

where labels e and o refer to the even and odd solu- 
tions respectively. In our model of ID dynamics we treat 
Eqs. (fTOf - fTTj) as boundary conditions for \z\ — * 0. 

In section llll Al we show that the phases tp e and ip are 
uniquely determined by l± and ip, and that their values 
can be calculated by solving the scattering problem in the 
quasi-lD geometry. This requires that l± 3> Ro, which 
allows to use the quantum-defect description of the short- 
range potential. As we discuss later, in contrast to 3D 
traps, in quasi- ID traps the values of ip are in general 
different for scattering waves of different symmetry. 



F. Time-dependent problem 

The atom-ion collision by moving the trapping po- 
tential d(t) is an intrinsically time-dependent problem 
which requires the integration of the time dependent 
Schrodinger equation for the given initial condition to 
predict the transition probabilities for the possible final 
states. In our approach to the dynamics we first calcu- 
late the correlation diagrams, showing the energy levels 
as a function of the trap separations, and on the ba- 
sis of these diagrams we predict the possible scenario of 
the atom-ion collision. When the motion of the trap is 
sufficiently adiabatic, the evolution of the system pro- 
ceeds along one of the energy curves, therefore the basis 
of adiabatic eigenstates, dependent parametrically on d, 
is particularly useful in the analysis of the collision pro- 
cess. Of course, the adiabaticity is usually broken in the 
vicinity of avoided crossings. In such cases, however, one 
can apply e.g. the Landau-Zener theory to calculate the 
probability of the adiabatic and diabatic passage through 
an avoided crossing. 

In our approach to the atom-ion collisions we consider 
that initially the traps are well separated. In this limit, 
the asymptotic states of ([T]) are given by products of 
the harmonic oscillator states in the two traps. In the 
course of dynamics, the distance d(t) decreases, particle 
interact for some definite time, and finally they are again 
separated, and the final state evolves into some superpo- 
sition of the harmonic oscillator states (see Fig. [1] for a 
schematic picture) . For such a scheme, we are interested 
in predicting the final state for some particular realiza- 
tion of d(t). We stress that in the intermediate phase, 
the system may evolve into atom-ion molecular complex, 
and such possibility is fully accounted for in our model. 

We start from the time-dependent Schrodinger equa- 
tion with the Hamiltonian |T|). We expand the time- 
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dependent wave function in the basis of the energy- 
ordered adiabatic eigenstates 4 , „(x 1 ,x 2 |d) of the Hamil- 
tonian (fT]) 

H{d)^ n ( Xi ,x a \d) = S n (d)* n ( Xi ,x a |d), (12) 

where we explicitly point out its dependence on the dis- 
tance d. Substituting the expansion 



*(Xi,X a ,t) 



2jc„(i) exp 

n 

x *„ (xi,x a |d(i)) 



-JdTE n {d{r)) 



(13) 



into the Schrodingcr equation, we obtain a set of coupled 
differential equations that govern the dynamics of the 
expansion coefficients c n : 



in = -d y^c m (^) exp 



dT (E n (d(T)) - E m (d(T))) 



X (*„(d)|^|*m(d)> 



(14) 



In the case of fully adiabatic evolution the coefficients 
c n (t) remain constant, and the evolution of the sys- 
tem proceeds along the adiabatic energy curves. From 
Eqs. (H3J) one can derive the condition for the adiabaticity 
of the transfer process. Adiabaticity requires that d mul- 
tiplied by the nonadiabatic coupling (^f n (d)\-^\^f m (d)} 
be much smaller than the frequency of the oscillating 
factor in the exponential of (Ti"4|) . In this case the oscillat- 
ing factor effectively cancels out the contribution due to 
changes of d. The adiabaticity condition can be written 
ashd(y„(d)\%\V m (d)) < (E n -E m ) 2 V m:7l 0. At trap 
separations where atom and ion can be approximately 
described by harmonic oscillator states, one can easily 
estimate the nonadiabatic couplings n (d)\^j\ 1 i' m (d)) , 
which gives the following constraint on the adiabatic 
changes of the trap separation: 



1<<{ E_ 
h h 2 u k 



E m y 



(15) 



where k stands for i (a) when the ion trap (atom trap) 
is moved. We stress that the latter condition is valid 
for transitions between different vibrational states, that 
may occur during the transfer of the atom or of the ion. 
In the case of avoided crossings between vibrational and 
molecular states, the adiabaticity of the transfer is char- 
acterized by the condition that can be determined from 
the Landau-Zener theory. 



III. RESULTS 

We turn now to the analysis of the adiabatic eigenen- 
ergies and eigenstates as a function of the trap separa- 
tion d, and of the trapping potential geometry. We start 
our analysis from the simplest case of quasi-lD traps, 
where the dynamics takes place effectively in ID, while 



the assumptions of equal trapping frequencies allows us 
to consider COM and relative motions separately. Before 
discussing this problem, we first study the dependence of 
the ID short-range phases ip e and ip on l±, R* and (p. 
We argue that in the considered range of parameters they 
are practically independent of the kinetic energy of the 
scattering particles, which is assumed in our quantum- 
defect-theory approach. The reader not interested in de- 
tails of the derivation may start at Section IIII B[ where 
we analyze ID relative motion eigenenergies and eigen- 
states assuming some particular values of tp e and tp Q . In 
Section fill CI we switch to 3D geometries, analyzing the 
system of two spherically symmetric traps with the same 
trapping frequency for atom and ion. Section [ill Dl is de- 
voted to the properties of adiabatic energy spectra in the 
vicinity of an avoided crossing. Using the semiclassical 
theory we calculate the level separation at the avoided 
crossing, and then applying Landau-Zener theory we in- 
vestigate the conditions for the adiabatic and diabatic 
transfer, depending on the distance between the traps. 
Finally, in Section Ull El we address the most complicated 
case of different trapping frequencies for atom and ion, 
coupling COM and relative degrees of freedom. Because 
of the complexity of this problem, we limit our analysis 
only to ID dynamics; we argue, however, that the ob- 
served behavior should be qualitatively valid also for the 
3D system. 



A. Short-range phases in quasi-lD traps 

As a preliminary technical step in preparation for the 
calculations of adiabatic energy curves that will be re- 
ported in the following, the present subsection deals with 
calculating the short-range phases ip e and ip in quasi- 
one-dimensional geometries. 

To find the ID short-range phases we solve the 
Schrodingcr equation for the relative part of the Hamil- 
tonian (fT]), assuming the same transverse trapping fre- 
quency for atom and ion: = luj_ ci = u>±, and ne- 
glecting axial trapping frequencies, which is valid for 
z <C Ri, R a - The latter condition requires Ri, R a ^> Rid, 
where Rid = max(R±,l±) determines the boundaries of 
the ID regime. In this way we obtain the following equa- 
tion for the relative wave function ^(r) 



h 2 „ 
2/i 



ae 
2^ 



1 2 2 



E *(r) = 0, (16) 



In the case of interaction potentials at distances Rid 
much larger than typical kinetic energies in the trapping 
potential: |V r (i?iD)| S> fiWi,huj a , we can neglect the en- 
ergy dependence of the short-range phases, and solve (|16p 
for E = fko±. 

At small distances, r <C Rid, the atom-ion interac- 
tion dominates over the transverse trapping potential, 
and the solution of (|16|) behaves according to At 
large distances, z 3> Rid, the motion in the transverse 
direction is frozen to its zero-point oscillations and in the 



8 



0.2 



0.0 



-0.2 



^0.05 
^ 0.00 



3D regime 

\ 

1 D regime 


• 


V(r)lh(o L 
K 

1D regime 


2 


1 

™ -8 








Y(z) ~ |z| sin(|z|/«*+p e ) 




1 0.09 /e* 



-1.0 



-0.5 



0.0 

zIR 



0.5 



1.0 



FIG. 4: (Color online) Potential for the atom-ion relative 
motion (upper panel) and the relative wave function (lower 
panel) for the quasi-one dimensional system with l± — 0.09R* 
and without trapping in the axial direction. 



asymptotic regime the scattering wave function assumes 
the form ^>(r) — > ?/>o(/ ) v f 'id(z) (\z\ — > oo), where ipo{p) 
is the ground-state wave function of the 2D harmonic os- 
cillator: -0o(/°) = e~ p /^'i)/^ 1 / 2 , and ^id(z) is a linear 
combination of odd and even waves: 

*l D (z) = C e ^ e (z) + C ^ (z), (17) 

with ^f e (z), ^g(z) given by ([TO]) and (jlip respectively 

Fig. 3] shows the potential for the relative motion, and 
the axial profile (p = 0) of the relative wave function for 
l± = 0.09i?*, which corresponds to a system of 40 Ca + 
and 87 Rb in a trap with ui± = 2n x 1MHz. In addition, 
the figure illustrates different length scales present in the 
quasi-lD problem, and indicates the ID (\z\ 3> -Rid) 
and 3D regimes (r <C -Rid) in the behavior of the wave 
function. 

Fig. [5] shows an example of the dependence of ip e and 
ip on (p, for ij_ = R* . It compares the results of nu- 
merical calculations with predictions based on the pseu- 
dopotcntial method, that is dicsussed in appendix IBl 
We observe that for even waves the agreement is fairly 
good, while for odd waves the agreement is poorer, which 
is probably due to the fact that the p-wave energy- 
dependent pseudopotential does not work already in the 
regime of R* ~ l± , or contributions of odd partial waves 
with / > 1 can be important. The pseudopotcntials are 
expected to give accurate predictions for R* <C l±- 

When R* 3> l±, the pseudopotential approach is not 
applicable at all, and one has to resort to numerical cal- 
culations. As an example, we present in Fig. [6] values 
for the short-range phases for l± = Q.1R*. We note the 
presence of several resonances in the dependence of tp e 
and <p Q . This behavior is related to the contribution of 
several partial waves for R* ^> lj_, leading to resonances 
when the energy of a bound state in the combined har- 
monic and r -4 potential becomes equal to the energy of 
the scattered wave. To analyze this issue more carefully 
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FIG. 5: (Color online) Even and odd short-range phases ip e 
and (fi calculated for l± = R* . Numerical results (red solid 
lines) are compared with predictions of the model replacing 
r~ 4 with the energy dependent pseudopotential. 



we have calculated the energies of some particular bound 
states in in the combined harmonic and r -4 potential. 
Fig.[7]shows the relation between the bound-state energy 
and the short-range phase for few bound states that are 
responsible for the scattering resonances close to <p = 
for (f> e (cf . Fig. [6|) . For comparison we include the bound 
states of pure r -4 potential and we label them by the 
angular momentum I. We note that, due to the limited 
resolution of our numerical calculations, Fig. [6] does not 
allow to resolve the single resonances close to = 0, 
however, one can observe the rapid changes of <p e in this 
region. 

For the numerical calculations we have transformed 
(|16[) into cylindrical coordinates and we were solving 
2D elliptic partial differential equation using a finite- 
element method. To fix the short-range phase we im- 
pose the boundary conditions on the logarithmic deriva- 
tive of *(r), 9 r *(r)/*(r), at r min , where r min = 0.09i?* 
(r min = 0.022i?*) for l ± = R* (l ± = 0.1R*) . At large 
distances we impose the Dirichlet boundary condition on 
the rectangle with boundaries \z\ = z max , p = p max , as- 
suming that the wave function has a Gaussian transverse 
profile at \z\ = z max and vanishes at p = p ma x- For 
l± = R* we took z max = 6.3i?*, p ma , x = 3i?*, while for 
l± = Q.1R* we have used z max = R* , p max = OAR*. 
To determine the values of tp e and </? , we fit at large 
distances (|z| 3> Rid) the symmetric and antisymmetric 
solutions of (fT6|) for E = fku± , to the asymptotic formula 

G3- 



B. Relative motion in a ID system: adiabatic 
eigenenergies and eigenstates 

In this section we consider the ID motion of an atom 
and an ion, assuming Ui = u a - In this case we can focus 
only on the relative motion described by the Hamilto- 
nian ([5]), with the boundary conditions at z — > 0, stated 
by (fTU|) and (fTTj) . To find eigenenergies and cigenfunc- 
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FIG. 6: (Color online) Even short-range phase <p e (upper 
panel) and odd short-range phase ip a (lower panel) calculated 
numerically for l± = 0.1R*, versus 3D short-range phase <p. 
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FIG. 7: (Color online) 3D short-range phase (p versus energy 
of bound states in the presence of the transverse confinement 
with l±_ = 0.1 J?*. The bound states in the combined harmonic 
and r~ 4 potentials (solid lines) are compared with the bound 
states of pure r -4 interaction (dashed lines). 



tions for arbitrary value of d we diagonalize the Hamil- 
tonian (|9|) numerically in the basis of its eigenstates for 
d = 0, that are found by numerical integration of the ID 
Schrodingcr equation. In the basis we include the low- 
est 50 odd and even eigenfunctions, which is sufficient 
to perform the diagonalization for d/R* < 2.5. Fig. [8] 
present the adiabatic energy spectrum for some example 
parameters: (p e = — 7r/4, ip a = ir/A and R* = 3.48/. The 
latter value corresponds to the system of 40 Ca + and 87 Rb 



in the trap with to = 2n x 100kHz, while the particular 
choice of the short-range phases (p e and ip Q is explained 
later in this section. Points with labels correspond to the 
wave functions shown in Fig. [5] presented together with 
their eigenenergies and potential energy curves. 

The three panels shown in Fig. [9] illustrate three dif- 
ferent regimes, where the system exhibits qualitatively 
different behavior. In the first regime, represented by 
the eigenstate ^ a and realized at large distances between 
traps, d 3> i? re i with R IC \ defined in (JSJ), the main effect of 
the atom-ion interaction is the distortion of the trapping 
potential, and in this limit the Hamiltonian can be diago- 
nalized analytically in the model of two coupled harmonic 
oscillators (see Section UlI El for more details). The value 
of the short-range phase is not important for this model. 
At large separations the eigenstate ^> a is typically only 
weakly perturbed with respect to the eigenstate of the 
harmonic oscillator. 

In the second regime, represented by the eigenstates 
and ^c, and realized at distances d ~ R Ie \, the system 
exhibits resonances between vibrational (\&i>) and molec- 
ular ( , F C ) states, manifesting themselves as avoided cross- 
ings in the adiabatic energy spectrum. The resonances 
appear when the energies of the two eigenstates become 
equal, and the energy splitting at the avoided crossing 
is proportional to the tunneling rate through the poten- 
tial barrier separating the two regions of the potential. 
In section UlIDI we calculate this splitting in the frame- 
work of a semiclassical approximation, and describe the 
dynamics in the vicinity of the avoided crossing. The 
value of the short-range phase determines the energy of 
the molecular states, but it is practically not important 
for the atom-ion vibrational states. Since the point of 
the avoided crossing depends on the short-range phases, 
the controlled collisions can provide some information on 
the short-range interaction potential. 

Finally, at distances d <C i? ro i the barrier separating 
the external trap from the well given by the atom-ion at- 
traction disappears and all eigenstates have some short- 
range component behaving at short distances according 
to Eqs. (fTUf and (fTTj) . Hence, in this regime all eigen- 
states depend on the value of the short-range phases. 
The molecular states with energies below the dissociation 
threshold (E = 0) are mainly localized in the well of the 
atom-ion attractive potential, however for the tight traps 
considered here, they can be strongly affected by the ex- 
ternal trapping (e.g. state , F e ). The other type of states, 
with energies E > 0, are analogs of the vibrational states 
and they are localized mainly at distances where the ex- 
ternal trapping potential dominates (e.g. state "F^)- 

So far we have discussed the properties of the energy 
spectrum for some particular choice of ip e and (p a . It 
turns out that a qualitatively similar behavior can be 
observed in all the systems with \ip e — ip Q \ = 7r/2. In 
the general case, however, the adiabatic energy spec- 
trum has a slightly more complicated structure, as it 
is illustrated in Fig. [TU1 showing the eigenenergies for 
(p e = ip = — 7r/4. We note that by going from d = 
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FIG. 8: (Color online) Energy spectrum of the relative motion 
for an atom and an ion in ID traps versus the distance d 
between traps, calculated for tp e — —n/4, tp a — 7r/4 and R* = 
3.48Z (see text for details). The dashed curve shows the line 
r;[iu! 2 d 2 , giving the approximate shift of the bound state in 
the trapping potential (see section MI CI for derivation) . 



to positive d ^ all the energies split into two branches. 
To understand the nature of this splitting in Fig. [IT] we 
present the wave functions for d = O.li?*, with the cor- 
responding eigenvalues marked with stars in Fig. [TDJ We 
observe that the branches represent the eigenstates local- 
ized on the left-hand side and right-hand side of the point 
3 = 0. The left-localized eigenstates correspond to the 
rising branches, because they are mainly localized in the 
region of strong atom-ion attraction. We note that this 
behavior is similar to the properties of eigenstates in a 
double-well potential, where appropriate wave functions 
are constructed by taking symmetric and antisymmetric 
combination of the states in two wells. In the particular 
case of tp e = ip , the symmetric and antisymmetric com- 
bination leads to the states localized on the positive and 
negative z semi axes, respectively. 



C. Relative motion in 3D for spherical u)t = u) a 
traps: adiabatic eigenenergies and eigenstates 

In this section we extend our analysis to 3D, and con- 
sider spherically symmetric traps for the atom and the 
ion with the same trapping frequencies: u>i = cu a . In 
this case the COM and relative motion can be separated, 
and in the following we focus only on the relative motion 
governed by the Hamiltonian ((4]). 

We have diagonalized the Hamiltonian (01 taking dif- 
ferent values of the short-range phase ip, and we have 
observed that the adiabatic energy spectra exhibit quali- 
tatively the same features as in ID. In the numerical cal- 
culations we first calculated the eigenstates for d = 0, by 
solving the radial Schrodingcr equation, for angular mo- 
menta I < 75 and for energies E < 75ftw. These states 
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FIG. 9: (Color online) Eigenstates of the relative motion for 
an atom and an ion in ID traps at different separations d 
between the traps, calculated for ip e = — 7r/4, ip = n/4 
and R* — 3.48/. Presented eigenstates correspond to labeled 
points in Fig. [8] The horizontal lines present the correspond- 
ing eigenenergies, and the thick black line shows the potential 
energy. 
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FIG. 10: (Color online) Energy spectrum of the relative mo- 
tion for trapped atom and ion in ID traps versus distance d 
between traps, calculated for ip e = cp a — — 7r/4 and R* = 3.48Z 
(see text for details). 



were used as a basis in the numerical diagonalization of 
the Hamiltonian (j4|). 

A sample adiabatic energy spectrum is shown in 
Fig. Q21 presenting the results for <p = — 7r/4 and R* = 
3.48/. At d = the angular momentum / is a good 
quantum number and the states have definite angular 
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FIG. 11: (Color online) Eigenstates of the relative motion for 
trapped atom and ion in ID at trap separation d — 0.1R* , 
calculated for ip e = tp = — 7r/4 and R* = 3.48/. The hori- 
zontal lines present the corresponding eigenenergies, and the 
thick black line shows the potential energy. 




FIG. 12: (Color online) Energy spectrum of the relative mo- 
tion for an atom and an ion confined in spherically symmet- 
ric traps versus the distance d between traps, calculated for 
Lp = — 7r/4, and R* = 3A8l r (see text for details). The blue 
dashed line shows the shift of the bound states with d, given 
approximately by ^/xu! 2 d 2 . 



symmetry, which is depicted by the appropriate symbols 
in Fig. 1121 For nonzero d, angular momentum is not 
conserved, hence in the dynamics the state can change 
its angular symmetry. We note that the energies of the 
bound states shift to a good approximation according to 
fiu> 2 d 2 /2. This behavior can be explained by noting that 
bound states l^moi) are concentrated around r = 0, and 
(*mol(d)|-ffrd(d)|*mol(d)) ~ #b + j^ 2 d 2 , where E b is 
the binding energy at d = 0. 



D. Avoided crossings: semiclassical analysis 

In the considered setup, level anticrossings reflect res- 
onances between molecular and trap states. A simple 
picture of such an avoided crossing is shown in Fig. 1131 
If one passes an avoided crossing from the direction of 
the asymptotic trap state, then for adiabatic evolution 
the system at small distances evolves into a molecular 
state. On the other hand, for diabatic passage, the 
particles remain in their traps, and the state basically 
does not change, apart from the modification due to the 
smaller trap separation. In addition, for fast changes of 
the trap positions, the particles can be excited to higher 
motional states. To describe quantitatively the dynamics 
in the vicinity of the avoided crossing one can apply the 
Landau-Zener theory. Assuming that close to the avoided 
crossing the eigenenergies are linear in d, and that d(t) 
varies linearly in time, the probability that the crossing 
is traversed diabatically is given by [lj, EB] 



P|i)—|i)' = CX P ~ 27r 



|(*i|ff|* 2 )| 2 
h\d\\dE 12 /dd\ 



(18) 



where the labels |1), |2) (|1)', |2)') respectively refer to 
the vibrational and molecular states before (after) pass- 
ing the avoided crossing (cf. Fig. [T3|) . and Ei 2 (d) = 
Ex{d) - E 2 {d). For K^liJltfa)! 2 « h\d\\dE 12 /dd\, 
the probability (fT8|) is close to unity, and the avoided 
crossing is passed diabatically. In the opposite case: 
|{*i|iJ|«- 2 )| 2 > H\d\\dE 12 /dd\, p is small, and the 
avoided crossing is traversed adiabatically. The matrix 
element (^i\H\^/ 2 ) can be related to the energy gap at 
the avoided crossing, and in this way we obtain the fol- 
lowing constraint on the adiabaticity of the transfer close 
to the avoided crossing: 



3E X dE 2 



dd dd 



< (£i(d) - E 2 {d)f 



(19) 



The position of the avoided crossings is directly related 
to the energies of the bound states; thus, a measurement 
of the final state after controlled collisions, provides tech- 
nique for spectroscopy of the trapped atom-ion complex. 
Since the avoided crossings become weaker as the separa- 
tion between traps increases, this scheme allows to probe 
only excited molecular states, having sufficiently small 
binding energies, comparable to the energy scales of the 
trapping potentials. 

In the case when the tunneling barrier is sufficiently 
large, the energy splitting at the avoided crossings can be 
estimated using a semiclassical approximation. We first 
focus on the relative motion in the ID system, assuming 
uji = cu a = lj. Using the standard WKB method, one can 
derive the following result for the energy splitting AE 

m 

AE = E^ moi (x 1 )^ vih (x2)V<x 1 )v(x 2 )Te- w / h . (20) 

Here X\, x 2 are some arbitrary points located in the 
classically forbidden region close to the classical turning 
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FIG. 13: (Color online) Schematic drawing of the avoided 
crossing between vibrational and molecular states in the sys- 
tem of a trapped atom and an ion. An adiabatic change of 
the distance d between traps induces a transition from the 
vibrational to the molecular states, while a diabatic process 
leaves the particles in the vibrational states. 



points, W(xi, X2,E) = dx^/2m{V(x) — E) is the ac- 
tion along the classical trajectory from x\ to cc 2 , T = 
is the tunneling time, 'Fmoi^i), ^vibO*^) denote, respec- 
tively, the molecular and vibrational states (cf. Fig. [T3|) 
with the same eigenvalues E (the resonance case), and 
v(x) = y / 2(V(x) — E)/m is the velocity of a particle with 
energy — E in the inverted potential — V(x). In princi- 
ple the choice of X\ and X2 is arbitrary, but in our cal- 
culations we take x\ and X2 at fixed distance from the 
location of the molecular (x = 0) and vibrational states 
(x = d). close to the classical turning points. In this 
way at sufficiently large separations ^moi^i), \E'vib(^2) 
becomes independent of d. 

We turn now to the case of relative motion in 3D traps. 
Similar to the analysis for ID systems, the energy split- 
ting at the avoided crossings can be calculated semiclassi- 
cally by applying the instanton technique. In our calcula- 
tion we adopt the formulation based on the path decom- 
position expansion developed by Auerbach and Kivelson 
[lq . We obtain the following formula describing the level 
splitting: 



AE = E^ mol (x 1 )* vib (x 2 )Vu(x 1 )?;(x 2 )Te- 



w/n 



(21) 



The meaning of all the quantities is the same as in the ID 
case, the only difference is the prcfactor A which accounts 
for the fluctuations around the instanton path (see (l6l | 
for details). In the case uji = uj a the instanton path that 
minimizes the classical action is simply a straight line 
connecting the trap centers. 

Fig. Q3] compares the semiclassical formulas Eqs. (|2"0)l 
and (f2~Tj) with the exact energy splitting determined from 
the numerical energy spectra, as those presented in Fig. [8] 
The presented results correspond to the avoided crossings 
between the vibrational ground state molecular states. In 
ID the energy splitting at the avoided crossing does not 
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FIG. 14: (Color online) Energy splitting at the avoided cross- 
ing between vibrational ground-state and the molecular state 
in ID and 3D systems. Semiclassical results (solid lines) are 
compared to the numerical values extracted from the adia- 
batic energy spectra for different combinations of ip e and tp - 
3D results present avoided crossings widths for s-wave bound 
states. 



depend on the symmetry of the molecular state, while 
in 3D it does depend on its angular momentum I, and 
we present here only the case 1 = 0. In the first ap- 
proximation, the splitting depends on the short-range 
phases only through the critical trap separation at which 
the resonance occurs. The numerical data are obtained 
for different combinations of the short-range phases, for 
which the adiabatic spectra exhibit avoided crossings at 
different values of d. In our approach, instead of calculat- 
ing ^moi(xi) and "Fvib^), we fix the overall amplitude 
in ([20|) by fitting to the single point at the largest trap 
separation, where we expect the WKB approximation to 
be most accurate. The semiclassical curves stop at the 
distances where the potential barrier disappears. We ob- 
serve that for the same ratio of R* to I the splittings in ID 
are larger than in 3D. Fig. [TH shows the prediction of the 
semiclassical formula Eq. ([20]) . and the energy splittings 
calculated numerically for molecular states with different 
I. We note that separations at the avoided crossing are 
largest for the spherically symmetric molecular states, 
with I = 0. 

The knowledge of the energy splitting AE can be used 
to calculate the probabilities of an adiabatic and dia- 
batic passage of avoided crossings. Assuming that the 
avoided crossing is traversed at constant rate, we apply 
the formula (H]) with dE x (d)/dd w and dE 2 (d)/dd = 
<*moi(d)l^rei(d)/ad|*moi(d)> « A^ 2 <2 (cf. Figs. 1 and 
[T2")) . which leads to 



P|i)-|i>' 



exp 



7T (AE) 2 

2h\d\nuj 2 d 



(22) 



where we use the same notation as in Fig. 1131 Analyzing 
Figs. [Inland 1 151 we can now estimate the rates d required 
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FIG. 15: (Color online) Energy splitting at the avoided cross- 
ing between vibrational ground state and molecular states 
with different angular momenta, for a 3D spherically sym- 
metric trap. Semiclassical results obtained by means of the 
instanton technique (solid lines) are compared to the numeri- 
cal values extracted from the adiabatic energy spectra. Semi- 
classical calculations stop at the distance d = 1.17?* where 
the potential barrier disappears. 



for adiabatic and diabatic transitions. For instance in 
3D, for the parameters of Fig. [T3J the diabatic transfer 
of particles across the avoided crossings up to distances 
dK, R* (AE < O.lhw) can be realized by keeping d/R* » 
O.OOlw . 

In summary, the analysis carried out here provides the 
basis for a description of the dynamics of various pro- 
cesses of interaction between an atom and an ion manip- 
ulated through external trapping potentials, and it gives 
a way to estimate quantitatively with simple analytical 
means the outcome of controlled interaction experiments 
in the different regimes described in Sect. II. 



E. Center of mass coupled to the relative motion: 
ID analysis for uii 7^ u a 

In this section we consider the effects of coupling be- 
tween COM and relative motions in ID, that appear 
when the trapping frequencies for atom and ion are not 
equal. 

Fig- HH shows the adiabatic levels as a function of trap 
separation for some example parameters: ioi = 5.5w a and 
l a = 0.97?*. This choice corresponds to the interaction 
of 40 Ca+ and 87 Rb in the traps with uj a = 2tt x 10kHz 
and LOi = 2n x 55kHz. To obtain the adiabatic spectrum 
presented in Fig. 1161 we performed the diagonalization 
of the Hamiltonian ^ in the product basis of the COM 
and relative motion eigenstates evaluated at d = 0. In 
the calculations we consider all the states with a total 
energy E < A60huj a , which leads to about 7100 states in 
the basis. 
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FIG. 16: Energy spectrum for atom and ion confined in har- 
monic traps with u>i = 5.5cj a , as a function of the distance 
d. Calculations are performed for tp e — —it/4:, ip = 7r/4 and 
l a = 0.97?* (see text for details). 



The arrows on the right-hand side of Fig. [TB] indi- 
cate the asymptotic states for large separations, which 
can be labeled by the number of phonons for the atom 
and ion trap, respectively. In comparison to the case 
where relative and COM motions are decoupled, we ob- 
serve the following new features: (i) the molecular spec- 
trum contains states with different numbers of excita- 
tions in the COM degree of freedom. This can be ob- 
served at d = 0, when the molecular levels are with a 
good approximation equally separated by Hu>cm where 
W CM = ( rn " UJ a + rrHuf)/M. [26|; (ii) the avoided cross- 
ings between molecular and vibrational states are weaker 
for the states involving vibrational excitations of the ion. 
This can be observed by comparing the avoided cross- 
ings for the state |0) n |l)i with the avoided crossings for 
the neighboring states. This behavior can be understood 
when we notice that for weaker atom trap (for oj a < oji), 
the particles have larger probability to tunnel when the 
atom is excited. 



F. Large distances between traps: two coupled 
oscillators approximation 

Finally we turn to the limit of large distances: d ^ 
Ri,R a . Expansion of the interaction up to second-order 
terms in the distance r leads to 



■ E 

u— i.a 



1 



2m,j 2 



l0—p-(z t - Zi)(z a - z a ), 



(23) 



where z v denotes the equilibrium position of the parti- 
cles and we neglect the second-order terms modifying the 
trapping frequencies. The Hamiltonian (|23p describes the 
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system of two coupled harmonic oscillators, which can be 
written in the form 

Hid = foi) a a<a + huj^b + hu c (a + a f )(6 + ft 1 ), (24) 

where we have introduced the usual annihilation and cre- 
ation operators 



2h 
2h 



Zi + i 



Pi 



Z a + I 



Pa 



m a uj a 



(25) 
(26) 



Here tu c denotes the coupling frequency: hio c = 
10E*(R*) 4 lJ t /d 6 . The validity of the considered model 
is limited by the assumption that a stable equilibrium 
position exists, which is fulfilled for Ri,R a < 0.57d, as 
it can be easily verified. Typical, maximal values of uj c 
for atoms in optical potentials and ions in rf traps are 
of the order of 10kHz. In practice the model ceases to 
be valid already at weaker conditions, when the terms 
higher than the second order cannot be neglected. In 
any case, this model is relevant for an important class of 
processes involving coherent transfer of quanta between 
the atom and the ion, of direct application in a quantum 
information processing context. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we analyzed in detail the interaction be- 
tween a single atom and a single ion guided by external 
trapping potentials. This kind of work is motivated by re- 
cently opened experimental possibilities within combined 
systems - currently being built in several groups world- 
wide - where magneto-optical traps or optical lattices 
for neutral atoms coexist with electromagnetic traps for 
ions. 

Tight confinement of single particles, associated with 
independent control of the respective confining poten- 
tials, allows for exploring different regimes of the two- 
body dynamics involving one atom and one ion at a time. 
At large distances, the interaction is dominated by the 
inverse quartic term arising from the attraction between 
the ion's charge and the electric dipole induced by it on 
the atomic electron wave function. A pseudopotential 
approximation, similar to that familiar from ultracold- 
atom collision theory, is not applicable in our case, as the 
characteristic range of the atom-ion molecular potential 
often exceeds the size of the tight traps experimentally 
available. To describe the interaction at short distances, 
smaller than that range, we employ quantum-defect the- 
ory, which allows to deal with different geometries from 
isotropic three-dimensional traps to very elongated quasi- 
onc-dimensional ones. 

A good description of adiabatic dynamics, involving 
processes where the traps are moved toward or across 
each other at a rate much slower than the trapping fre- 
quencies, can be obtained from quasistatic cigenencrgy 



curves, which we calculate for various trapping configu- 
rations based on the methods outlined above. A remark- 
able feature displayed by the system is the presence of 
resonances between molecular-ion bound states and mo- 
tional excitations within the trap. These trap-induced 
resonances are similar in nature to Feshbach resonances 
driven by external fields, and they could as well be spec- 
troscopically probed in simple experiments where the in- 
teraction is controlled via the external guiding potentials. 

In addition to the above aspects, an important moti- 
vation for the interest in systems of trapped atoms and 
ions and their trap- induced resonances resides in possible 
applications to quantum information processing. In this 
context, one can utilize controlled atom-ion interactions 
to effect coherent transfer of qubits, thereby creating in- 
terfaces between atoms and ions. By storing quantum in- 
formation in internal atomic states, and performing gate 
operations with ions, one would combine the advantages 
of both: (i) long decoherence times for neutral atoms 
and (ii) short gate-operation times for charged particles 
due to the relatively strong interactions. In the present 
work we have focused on the motional degrees of free- 
dom, which can serve as auxiliary degrees of freedom for 
quantum gates involving internal-state qubits, provided 
an appropriate coupling mechanism between internal and 
external degrees of freedom is employed - e.g., sideband 
excitation via a laser. 

Another potential application of our results is cooling 
of the atomic motion. In typical neutral-atom quantum 
computation schemes with qubits stored in internal states 
of atoms trapped in optical-lattice sites, motional excita- 
tions constitute a serious source of errors, and the atomic 
motion needs to be cooled in a state-insensitive manner 
between computational steps in order to avoid qubit de- 
coherence. Since the long-range part of the atom-ion in- 
teraction is not sensitive to the internal state [27| . our 
setup can be applied for sympathetic cooling of atoms, 
through the exchange of energy with laser cooled ions 

Beside these applications in the context of quantum 
information processing, our work opens broader perspec- 
tives for the study of new interesting collisional physics 
in a physical situation never explored before. In prin- 
ciple, our scheme allows for production of ultracold 
trapped atom- ion molecules (molecular complexes), when 
the trapping potentials are lowered adiabatically at the 
stage when the particles remain close to each other. In 
this way our method can be regarded as a way to per- 
form cold chemical reactions, where the final state of 
the molecule can be well controlled. Beyond sufficiently 
large separation between traps, however, the survival 
probability of such molecules in the final state is neg- 
ligible. Indeed, in the calculations presented here, we 
have only regarded processes leading to a final state in 
which the atom and the ion are still separated, and no 
molecule has been formed as an outcome of the inter- 
action. In other words, the Hamiltonian describing our 
single-channel model docs not include the possibility of 
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transitions to other Born-Oppcnheimer curves. While 
this is a justified assumption for the dynamics considered 
here, it is in a sense a limitation of our current approach. 
Investigation of molecular formation in traps beyond this 
approximation is certainly among the interesting devel- 
opments that can arise from our work. Its long-term 
motivation is really to open, beyond the present exam- 
ples, a new paradigm for cold collision physics, which can 
be described as the mechanical control of single-particle 
chemical reactions. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN Q 

While the Hamiltonian ((T|) is intuitively obvious, we 
find it nonetheless worthwhile to summarize the micro- 
scopic derivation in an adiabatic approximation, and to 
discuss its validity. For simplicity we consider an atom 
with a single outer-shell electron (alkali atom) and an 
ion with single positive charge +e. In addition we do 
not consider the internal structure of the ion. The total 
Hamiltonian can be written as 



H 



PI 



rr 
P2 



2mi 2m 2 



H c + H } 



las 



(Al) 



Here the labels 1 and 2 correspond, respectively, to 
atomic nucleus and ion, and H c is the Hamiltonian of 
the electron, which includes the Coulomb interactions. 



2m c 



x e - Xl 



X e - X 2 



X 2 - Xl 



(A2) 



Here, for simplicity, we omit the contributions of the 
core regions for both the atom and ion. In this way our 
model refers in fact to Hj molecule. For alkali atoms and 
alkaline-earth ions one should treat the complete struc- 
ture of core regions, however, our approach can be readily 
generalized to this more complicated case. 

The term Hi as describes the interaction of the atom 
with a laser beam creating an optical potential, written 
in the electric dipolc representation 



where d = e(xi — x e ) is the dipole moment of the atom, 
and for simplicity we neglected the influence of the laser 
on the ion, which has typically a different electronic 
structure that the atom. In addition we have applied 
the long-wavelength approximation, neglecting changes 
of the electric field on the scale of the atom. In the case 
of the optical lattice the transverse part of the electric 
field E_l can be assumed to have the form of a standing 



Ej_(x, t) = E cos(ui L t) [cos(k L x + (j) x ) + (x -> y, z)] , 

_(A4) 

where for simplicity we have assumed the same amplitude 
Eo, the same wave vector fc^, and the same frequency of 
the laser light uj, for all three laser beams creating the 
optical lattice potential. The abbreviation (x — ► y, z) 
denotes sum of the terms with x replaced by y and z, 
and (fk for fc = x, y, z is the phase factor characterizing 
the standing wave. Finally, H r { is the electric potential 
creating the rf-trap 



H Ti = e$(xi,£) + e$(x 2 ,i) - e$(x e ,t), 



(A5) 



where $(xi,i) is the time-dependent electric field of the 
rf trap 

$(x, t) =^(u x x 2 + u y y 2 + u z z 2 ) 



+ —{v x x 2 + v y y 2 + v z z 2 ) cosa; r ft. (A6) 



-dE ± ( Xl ,t), 



(A3) 



Here, w r f is the frequency of the time-dependent part of 
the electric potential, and Uk, Vk (k = x,y,z) are am- 
plitudes depending on the geometry of the trap [f| ■ The 
electric field at every instant of time has to fulfill the 
Laplace equation: A$ = 0, hence, the coefficients Uk, «fe 
are subject to the following conditions: u x + u y + u z = 0, 

V X + Vy + V z = 0. 

Below we indicate the basic steps of the derivation. 

Expansion in the basis of Born-Oppenheimer wave 
functions for electron motion. We start from generating 
a complete set of electronic wave functions, parametrized 
by the positions of the atomic core and of the ion 

fr e $ n (x e |xi,x 2 ) = E n (\x 2 -xi|)$„(x e |xi,x 2 ) (A7) 

In this way the total wave function can be expanded in 
the basis of Born-Oppcnheimer electronic wave functions 

*(xi,x 2 ,x e ,t) = ^c n (xi,x 2 ,i)$ n (x e |xi,x 2 ). (A8) 



Since the basis is complete, the expansion of the wave 
function does not involve any approximations. 

Retaining in the expansion only the modes coupled by 
the laser. In the expansion (|A8|) we keep only two modes 
coupled by the laser light creating the optical lattice: the 
electronic ground state $ s (x e |xi,X2) and the electronic 
excited state 3> e (x e |xi, x 2 ). In this way we utilize the 
Born-Oppcnheimer approximation, treating the electron 
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motion in the adiabatic approximation. This assumes 
that the time scale of the electron dynamics is much 
faster than the dynamics of the atomic nucleus and of the 
ion, which is typically fulfilled since the electron is much 
lighter than the other two particles. The approximation 
of the two coupled channels can be easily generalized to 
more, or even infinite number of channels, since the other 



channels arc weakly populated and we treat them within 
the perturbation theory. 

Adiabatic elimination of the excited electronic state 
coupled through the laser: derivation of the optical trap 
potential. The expansion coefficients c 5 (xi,x 2 ,t) and 
c e (xi,X2,t) fulfill the following set of coupled equations: 



ih 
ih 



d c g (xi,x 2 ,t) _ f _p|_ 
dt 

5c e (xi,x 2 ,t) 



\2mi 2m 2 



E g (\x2 - Xi|) + e<I>(x 2 ,i) c s (xi,x 2 , t) - d C gEj.(xi, i)c e (xi, x 2 , i) 



dt 



P? 



rr 
P2 



2toi 2m 2 



(A9) 



E e (\x 2 -xi|) + e$(x 2 ,t) c e (xi,x 2 ,t) -degEj_(xi,t)c 9 (x 1) x 2 ,t), (A10) 



where d C g(xi,x 2 ) = ($ e |d|<I> s } is the dipole matrix ele- 
ment between the ground and excited electronic states, 
which in general depends on the position of atom and 
ion. In the derivation of (|A9p - (|A10p we have neglected 
the action of the rf field on the atom core and the elec- 
tron, putting 



(*fc|fr rf |*fc) 
($ e |ff rf |* fl ) 



e$(x 2 ,i), k 
0. 



(All) 
(A12) 



We focus on the regime of far-detuned laser: A ^> fJ^, 
where A denotes the detuning A(|x 2 — xi|) = i? e (|x 2 — 
xi|) — E g (\x2 — xi|) — Ul, and = — d cg E /Zi is the the 
Rabi frequency. For such conditions, the excited state 
is only weakly populated and can be adiabatically elimi- 
nated. Additional simplification comes from the fact that 
the transitions between states g and e, due to the laser 
light, occur on a time scale much shorter than the mo- 
tion of atom and ion, and the dynamics of the atom and 
ion centers of masses can be decoupled from the internal 
dynamics, in full analogy to the Born-Oppcnhcimer ap- 
proximation. Eliminating the excited state, in basically 
the same manner as in the standard derivation of the AC 
Stark shift, we obtain the following equation that governs 
the dynamics of the atom and ion center of masses: 



be justified, since in the atomic collisions, the probabil- 
ity of finding the particles at short distances, comparable 
to the range of chemical binding forces, is typically very 
small. 

Time averaging over fast oscillations of the rf field: 
derivation of the Paul trapping potential. We replace the 
time-dependent rf field by an effective, adiabatic poten- 
tial. In this way we neglect the fast micromotion of the 
ion on the time scale of w r f 01 • Finally we obtain the fol- 
lowing equation that describes the dynamics of the wave 
function c g (xi,x 2 ,i) dependent on the atom and ion po- 
sitions 



dt 



PI 



r> 2 

P2 



2mi 2to 2 



£ ff (|X2-Xi|) + Vopt(xi) 



y r f(X 2 ) Cg(Xl,X 2 ,t). 



Here V T f denotes the effective potential 

Kf (x) = -m 2 (uj x x 2 + ujyy 2 + uj z z 2 ) 



(A15) 



(A16) 



ih 



<9c s (xi,x 2 ,i) _ / p 2 



r> 2 

P2 



dt 



\2toi 2m 2 



E g (\x 2 -x 1 |) + y opt (x 1 ) 



where 
V r op t(x) 



m 2 L 
1a 



e$(x 2 ,t) )c g (xi,x 2 ,t), (A13) 



(cos(fc L x + 4> x f + (x y, zj) (A14) 



is the effective potential due to the laser field. In the for- 
mula (|AT4)) we have neglected the position dependence 
of the detuning A(|x 2 — xjj) and of the dipole matrix el- 
ement d cg (x 1 ,x 2 ), assuming that they are modified only 
at very short distances between the particles. This can 



where uj k = (a k + ^/2) 1 / 2 w rf /2, a k = 4eu fe /(m 2 w 2 f ) 
and qk = 2ew/ i: /(m 2 w 2 f ) for k = x,y,z. The interac- 
tion between atom and ion is given by i? s (|x 2 — xi|) - 
the ground-state energy of the electron motion, which at 
large distances behaves as E g (r) ~ — ae 2 /(2r 4 ). 



APPENDIX B: DETERMINING OF <p e AND p 
FROM PSEUDOPOTENTIALS 



In principle, the replacement of the r~ 4 interaction 
by the pseudopotential is strictly valid when R* <C lx- 
Nevertheless, it is possible to try to determine ip e and ip a 
for R* ~ Z* using the energy-dependent pseudopotentials. 
To this end we replace the r -4 interaction with [l7|, IH| 



2irh 2 a(k) d 

V s (r) = <Hr)^-r 

/i dr 



(Bl) 
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for even scattering, and [lj 



V P (r) 



nh 2 a p (k) s 



(B2) 



In the pseudopotcntial approximation one can solve the 
quasi-lD scattering problem exactly, and calculate values 
of the ID scattering lengths for even [13, [U and odd [HJ 

waves 



for odd scattering. Here, a(k) = — tern So(k) / k is the 
energy-dependent s-wave scattering length, a p is the p- 
wave scattering length: a p (k) 3 = — tan<5i(fc)/fc 3 , the 

symbol V (V) denotes the gradient operator that acts to 
the left (right) of the pseudopotential, and So(k), 6%(k) 
are the s- and p-wave phase shifts, respectively. 

In ID the even (ip e ) and odd (i/j ) scattering waves have 
the following asymptotic behavior: ipe(z, k) ~ sin(fc|z + 
S e (k)), and ip (z,k) ~ z/\z\ sin(fc|z + S Q (k)) (\z\ — > oo), 
where 5 e (k) and 8 (k) denote even and odd scattering 
phases, respectively. In analogy to the three-dimensional 
scattering theory, we define ID even (af D ) and ID odd 
(a° D ) scattering lengths: = limfc_>o — (tanS e ,o(k))/k. 
From this definition it follows that for k = the even and 
odd scattering waves behaves as: ipe(z, k = 0) ~ \z 



and i/) (z,k — 0) ~ z 
result applied to Eqs. 



1D z/|z| (\z\ -> oo). 
and (fTTj) leads to 



'ID 



a° m /R* 



/R* = — COt ip e 



COt (fi 



u, 1D 

The latter 

(B3) 
(B4) 



'ID 



(fc) = 



2a s {E) 



1 3 

2'2 



«!d(*0 = f 



12a p (E)'' 



_ 1 3__E_ 

2' 2 2fej ± 



(B5) 



(B6) 



where = huj± + h 2 k 2 /(2/x) and C(s, a) denotes the Hur- 
witz Zeta function: C( s i a ) = S^Lo(^' + a ) _S Fi- 
nally, to relate the energy-dependent scattering lenghts 
a s (E) and a p (E) to ip, we apply the exact solutions of 
the Schrodinger equation for r -4 potential, given by the 
Mathieu functions ^BQ- 
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